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Abstract 

We study relaxational behavior from a highly excited state for a composite system 
in partial contact with a heat bath, motivated by an experimental report of long-term 
energy storage in protein molecules. The system consists of two coupled elements: The 
first element is in direct contact with a heat bath, while the second element interacts only 
with the first element. Due to this indirect contact with the heat bath, energy injected into 
the second element dissipates very slowly, according to a power law, whereas that injected 
into the first one exhibits exponential dissipation. The relaxation equation describing this 
dissipation is obtained analytically for both the underdamped and overdamped limits. 
Numerical confirmation is given for both cases. 

1 Introduction 

When a system is in contact with a heat bath, in the most common situation it relaxes toward 
the equilibrium state rather rapidly, typically in a manner characterized by exponential decay, 
except at a phase transition point. For instance, a hot particle in cold water generally loses its 
energy very quickly. However, does such rapid relaxation of a high-energy state hold generally? 
Consider, for example, a system composed of several elements that is in contact with a heat 
bath only through its outer part. A protein molecule in water represents such an example, 
since the hydrophobic part of a protein molecule will be segregated from an aqueous solution 
[|TJ. In such a situation, the relaxation process of the internal part naturally depends on the 
dynamics of the whole system, including the temporal and spatial relationship between the 
internal and external parts. Indeed, it has been reported that proteins excited into high-energy 
states, exhibit very slow relaxation, compared with that seen in standard dissipation processes 
in systems of similar size 0. Such long-term energy storage may be necessary for proteins to 
work functionary. For instance, consider a thermal ratchet model [§, f|, || for motor proteins 
that work in a heat bath. In the model, the existence of both high and low temperature domains 
in proteins are assumed, to generate a directive motion. Then, for such a mechanism to work, 
it should be answered how the proteins maintain a high energy (i.e., keep an effectively high 
temperature state) at a localized domain within them. 
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In connection with this long-term energy storage in proteins, two of the present authors 
PI have studied a model of a coupled pendulum in partial contact with a heat bath. It was 
found that when a pendulum that is not in contact with a heat bath is highly excited, the 
energy dissipates only slowly, approximately with logarithmic decay. This slow relaxation 
seems to be consistent with the Boltzmann- Jeans conjecture J7|, ||| for very slow relaxation in 
Hamiltonian systems without contact to a heat bath. In such systems, the relaxation time 
increases exponentially of some positive power of the injected energy || |J. 

When all elements in a Hamiltonian system are in full contact with a heat bath, the mech- 
anism allowing for slow relaxation may be destroyed completely. On the other hand, when a 
system contacts a heat bath only through an external part, the internal part can evolve ac- 
cording to its 'own' dynamics. For this reason, the internal Hamiltonian dynamics may have a 
strong influence on the relaxation process, and this can lead to a slow relaxation process and 
long-term energy storage. It is thus important to study how a system in partial contact with a 
heat bath insulates its internal part from the effect of the heat bath. 

To study such behavior of Hamiltonian systems in partial contact with a heat bath, we 
adopt a simple model, introduced in §2, consisting of two coupled elements: The first element 
is in direct contacted with a heat bath to allow for the dissipation of the energy, while the 
second element interacts only with the first element, and not directly with the heat bath. 
We study relaxation from a far-from-equilibrium state by preparing the second element in a 
highly excited state with a large kinetic energy. The first element, on the other hand, remains 
always in near-equilibrium with the heat bath. It should be noted that, without the heat bath, 
the presently studied Hamiltonian system is always integrable and does not relax, even when 
prepared in a highly excited state. Therefore the relaxation of the second element to thermal 
equilibrium is due to the indirect effect of the heat bath. In §3, the relaxation equation for the 
underdamped limit is obtained, while the case of the overdamped limit is studied in §4. The 
validity of these equations is confirmed by numerical simulations discussed in §5. Relaxation 
in the intermediate damping case is examined in §6. The results of numerical experiments for 
this case are also reported there. From these analytical and computational studies, two kinds 
of power law are obtained for the underdamped and overdamped cases, respectively. It is found 
that these power laws are independent of the temperature. 



2 Model System 



To study relaxation of a system with internal degrees of freedom that is in partial contact with 
a heat bath, we consider a Hamiltonian system with 2 degrees of freedom of the general form 



H 



(1) 



where U{6) is a periodic function of period 2n satisfying \U\ < oo. Considering the system to 
be in partial contact with a heat bath, and adopting the standard representation of the heat 
bath, we take the motion of the elements to be given by the following Langevin equations: 

9U 

Pi = ~— ~7Pi + £(*), 

(2) 



P2 



dU 

Wo 
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Here the random force £(t) due to the heat bath consists of a Gaussian white noise satisfying 
= an d (^(^1)^(^2)) = 2jT5(t2 — tx), with (•) representing the temporal average. (We 
use units in which ks = 1). The parameter 7 is the dissipation coefficient characterizing the 
dissipation to the heat bath, and T is the temperature of the heat bath. With such contact to 
a heat bath, it is generally expected that a system of this kind will relax to equilibrium. The 
change of the total energy E due to dissipation to the heat bath can be written 

E = - 1P l+p 1 Z(t), (3) 

and we see that E remains constant if there is no contact (7 = £ = 0) with the heat bath. 

Using the set of variables = Q\ + 9 2 , <fi = #1 _ #2 an d P = px + P2 = O, we rewrite the 
Hamiltonian as 

P 2 6 2 

H = — + ^ + U{<P). (4) 

When there is no heat bath (i.e., for 7 = £ = 0), the equation of motion is completely separable 
in these variables: P = 0, <fi = —2dU/d<p. Thus we see that in this case this system is completely 
integrable, and the total momentum P is a conserved quantity. 

In the following sections, we study this system analytically to elucidate the relaxational 
behavior in certain limits. The results of numerical simulations of some models are given in §5 
and 86. 



3 Underdamped Limit 

In this section we consider (2) in the underdamped limit, i.e. with a small value of 7. With 
the presence of this weak damping term, the total energy E and total momentum P, which are 
conserved quantities for 7 = 0, become slow variables. 

3.1 Derivation of relaxation equation 

Let us consider the situation in which the system is initially far from equilibrium, with very large 
IP2I 3> \pi\- The rotation of 0, accordingly, is much faster than its change due to dissipation. 
Then, the terms dU/ d9x and dU / d8 2 in (|J) change so fast that their effects on the dynamics of 
Px and P2 are averaged out over a time scale relevant for these slower variables. Here we have to 
distinguish the orders of characteristic time scales. We have the following relations : (the time 
scale of £(£)) <C (the time scale of rotation 27r|</>|~ 1 ) (the time scale of relaxation 7 _1 ). Now, 
we introduce a slower time scale r satisfying 27r|</>|~ 1 <ti r 7 _1 , and consider averages over 
this time scale, X(t) = - // +r Xdt. Using variables averaged over this time scale, the mode of 
fast rotation is eliminated. 

The equations of motion for the slow time scale are obtained by integrating (Q) from t to 
t + r. This yields 

t ^ "--^ 

(5) 



dp2 
dt 



U> 
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where U' = df//d0 = 0U/d6 1 = -dU/d6 2 , U' = \ J t t+T d[/(0(£'))/d0dt' and (f(*)£(f)) = 
2^T5(t — t'), with = 0. The 5-function here is defined with respect to the time scale r. 

In the limit \<f)\ — > oo, the averaged variables can simply be replaced by their uniform 
averages over 0, defined generally by X = J 27r X(0)d0/(27r). Then we obtain 



dpi _ T 



(6) 



where U' — for periodic potentials. Then, as the lowest-order estimation in this limit, we 
have 

Pi-p^e-t 1 , j5 2 =p 2 (0), (7) 

starting in a state far from equilibrium with large values of \pi\ and |p 2 |. Beginning in such a 
state, after a time scale of 0(7 _1 ), the first element approaches thermal equilibrium (pi ~ T), 
while, within this lowest-order approximation, the second element continues to satisfy p 2 = £>2(0) 
and does not evolve toward equilibrium. 

To observe the relaxation process of the second element, we need the next order estimation. 
Since is finite and can depend on 0, in general the value of a time-averaged variable X deviates 
from that of the corresponding phase-averaged variable X. We now derive the correction to the 
equations of motion that this deviation introduces to first order in \<fi\~ l , and thereby obtain 
the relaxational equations for the second element and the total energy E. 

After the exponential decay derived by (^) occurring on a time scale 0(7 _1 ), p\ is approxi- 
mately equal to its equilibrium value while p 2 is still large. Then the relations 

|pi| ~ Vr<^p 2 ~P (8) 

hold, where P is the total momentum of the system and is large enough for the fast rotational 
motion of the second element to be maintained. From (|8]), we straightforwardly obtain <fi ~ —P. 
The total energy E of the system can then be expressed by the variables P, <p and p\ as 

E = ^- + U( ( f>)-p 1 P + pl (9) 

where U(<f)), p\P and p\ are much smaller than P 2 /2. 

Since |0| is sufficiently large and 7 is sufficiently small, E and P are almost conserved over 
the time scale r. To consider the slow evolution of E and P, we eliminate the fast variable 0, 
using the time averaged variables X. Within the slow time scale, the value of E can then be 
expressed as 

E = ^- + U-p 1 P+pl (10) 

where - f t +T p\dt' is approximated by p\ on the time scale r 7 -1 . From the relation between 
(f|) and ([K]), the value of is estimated as 




-), (11) 
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where 5U = U — U . The last expansion is valid if the conditions \SU\ <C P 2 and T ^ P 2 are 
satisfied. The momenta of the two elements are then given by 

P + 5U 2p 1 5U ^ f 5U 2 T5U . , , 

This result reveals that both elements oscillate rapidly around their averages pi and £>2 (= 
P — pi). We see that the amplitude of oscillation becomes smaller as the second element 
rotates faster (i.e., for larger p 2 and P). 

In order to obtain the relaxation equation for P, we substitute (12 ) into By considering 



only the case with a sufficiently low temperature, and neglecting terms of 0{TP ), we obtain 
• 78U 2 2p 1 5U 6U£ „. „ ^ ^,5U 2 T . 

^ = - V" + p + p + Pl(_7Pl + + 0{ ~W> pi y (13) 

This equation includes fast variables, while the total energy P is a slow variable. It is thus 
appropriate to average out the fast change in this equation. Doing so, we obtain 

dP t 1 r t+T n 11 r t+T d v 2 ~ 5U 2 T 

X—fr-rl SUdt '^P^l WO + ^I + O + O^,^), (14) 



by using (||). To first order in \<p\ 1 , we can replace the time average of 5U 2 by the cycle average 

- ( t+T 5U(<p(t')) 2 dt' -> — r 5U{<p) 2 d<j> = W 2 . 
t Jt 2ix Jo 



5U 2 : 



Applying this approximation, we obtain 

dP ^W 2 1 r t+T CTT ^ nil d,p? ~. ^ r 5U 2 T . 

Since E = U + pf/2 + p\/2 and \p%\ ^> \pi\ (~ vT), ( |15l) can be rewritten in terms of 
^4 i ,.*4-<r xrr2 



-(^f) = -rySlP+M- / ^')dt'} + 0(^-,T). (16) 
dt 4 T Jt p 2 



This equation describes the relaxation to thermal equilibrium for the second element. To 
consider only the long-time (3> 0(t)) behavior at the current order of approximation, it can 
be written as 

±(£) = ^§U* + p 2 6UZ(t). (17) 
3.2 Relaxation at T=0 

When T = 0, £(i) vanishes, and pi is negligible (or, more precisely speaking, it is 0(P -3 )). 
Then we obtain the following relations from flTTD and (O) 

• _ 2SU 5U _ 5C/ 

= -P + -p-, Pi = -p, P2 = P~y. 



The relaxation equation for the total energy P, (|15|), is then simplified as 

dP -fWP 



it 2E ' < 18 > 



for U <C E (that is E ~ P 2 /2). The relaxation equation for the second element is now given 
by 

^(f) = -7^- (19) 

These relaxation equations reveal that a coupled system in partial contact with a heat bath 
exhibits non-exponential, slow relaxation after excitation to a high-energy state. Indeed, from 
flTSP and (0), the form of relaxation is obtained as 



E = E^-iSlPt, (20) 
p 2 (tf = p 2 (0) 4 -4 7 ^t, (21) 



where E and P2{0) are the initial total energy and momentum of the second element at t = 0, 
which is defined as a time at which the first element reaches thermal equilibrium. Thus, the 
characteristic time t r for the relaxation of p 2 is given by 

P2(0) 4 El . . 

With an increase in the injected energy (or momentum), the relaxation time increases in propor- 
tion to £>2(0) 4 . This increase is slower than the exp(p 2 (0)) form obtained from the Boltzmann- 
Jeans conjecture ||, but even here, there is a rapid increase of t r with the injected energy. 
Comparing the characteristic time (|22"1) to that of the standard dissipation processes exhibited 
by systems in full contact with a heat bath (where t r ~ 7 _1 ), it is confirmed that partial contact 
indeed results in long-time relaxation behavior. 

3.3 Relaxation at finite temperature 

Using the expression (|TTD, we derive a relaxation equation for finite temperature. The solution 
of fll7D can be written 



P2{tf P2W 



-j5U 2 t+ C p 2 5U£(t')dt' (23) 
Jo 



for t 3> 2n\p2\~ 1 . We next take an ensemble average (indicated by (•)) of (p3|), in order to 
see the averaged behavior for the long-time profile of the relaxation. Although the ensemble 
average of the third term, (p 2 SU^), cannot necessarily be decomposed into the product of the 
average of each quantity, such a decomposition is valid at the current order of approximation, 
as discussed in Appendix 1. Hence we get 

m« K Mm ^ lW -n + jy V)mmmV )) A i = <pfi _ lWPt , (24) 

Since E ~ p\j1 for \p 2 \ 3> this result leads to 



(E 2 ) = (E'}- 1 6Un. (25) 

Thus the relaxations of the mean values of E 2 and p\ exhibit the same forms as these in the case 
with T = 0. The validity of the approximations used to obtain these forms have been confirmed 



numerically, as discussed in §5, where ( |25| ) to be valid over a sufficiently long interval. 
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Furthermore, we can approximate the evolution of the variance of E 2 (or p|) around the 
mean value, 

((E 2 - (E 2 )) 2 ) ~ f / < (p 2 (t')P2(t")^(t , )^(i"))^(i , )e(t"))dt'dt" = 2 7 T f \ P2 {t') 2 5U{t') 2 )dt' , 

JO J JO 

(26) 

where we have again used the decomposition discussed in Appendix 1. If the time t is not too 
large, this variance can be approximated by 

{(E 2 -(E 2 )) 2 } ~2 1 T(p 2 (0) 2 }M 2 t~4 1 T(E )M 2 t. (27) 

This form for the increase of the variance has also be shown to agree with numerical results, as 
discussed in §5. 

From (25), we obtain the mean relaxation time (t r ) as 

< ( -> - H- (28) 

We thus find that the value of the temperature has little effect on the nature of the relaxation 
from a high-energy state, as long as the temperature is sufficiently small. 



4 Overdamped Limit 

When the dissipation is sufficiently large, the inertial force can be neglected. In this overdamped 
limit, the equations of motion (Q) can be written 

8TI 

pi = -r l Wi +r l m+o^- 2 ), 

(29) 

dU 

P2 = ~df 2 

Then, for sufficiently large 7, the equation for the energy dissipation is obtained straightfor- 
wardly as 

E = - 1V \ + Pl £(i) = -r'U' 2 + 7- V£(i) + 0( 7 - 2 ), (30) 
where U' = dll/d(j) = dU /dd\ = —dU /d9 2 - Solving this equation, we have 

E = Eq — 7- 1 f U' 2 dt' + T 1 f U'^t')dt', (31) 
Jo jo 

with the initial energy Eq at t = 0. 

We now consider the case 7 3> P/(2n). In other words, that in which the time scale of 
dissipation from the first element is much faster than the rapid rotation of the second element. 
Here, we are interested in relaxation from a high-energy state (with very large \p 2 \) taking place 
over a much longer time scale than the rotation time (2irP~ 1 ). In this case, the integration over 
time can be approximated by the cycle average: | J * U' 2 dt' ~ J 27r U' 2 d<f) = U' 2 . By using this 
approximation, (|3l|) can be replaced by 

E = E - 7 - 1 C7^i + 7' 1 /' U'i{t')dt'. (32) 

jo 

An intuitive explanation for this approximation is as follows: In the limit P — > 00, the change 
of 4> is expected to be uniform in time. Then relaxation on the slower time scale can be 
approximated by using the uniform average of U'(<f)) over 0, to lowest order in 7 -1 and P^ 1 . 
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4.1 Relaxation at T=0 



For T = 0, there is no random force, and the relaxation is obtained directly as 

E = E - 7- 1 £7 75 t. (33) 

Since E ~ p|/2, this leads to relaxation of the second element given by 

p 2 (t) 2 =p 2 (0) 2 -2 1 - 1 U^t, (34) 

in which P2{0) is the initial momentum. 

From ([34]), the characteristic time t r for the relaxation is obtained as 

t r = ™!E = 25. (35) 

Thus in the present case, the dependence of t r on p2{0) is different from that in the underdamped 
case, where t r ~ p 2 (0) 4 . We remark that even in this overdamped case, the relaxation is much 
slower than in the case of full contact with the heat bath (t r ~ 7 -1 )- The crossover between 
the two limiting cases is discussed on the basis of numerical computations in §6. 

4.2 Relaxation at finite temperature 

Taking an ensemble average of fl52]), we consider the relaxation of the mean value of E: 

(E) ~ (E ) - r X U^t + I' 1 l\u')(mW = (Eo) - r'U^t. (36) 



Here ({7'(</>)£(£')) has been decomposed as (U'(<fi))(£,(t')), as is justified at the present order 
of approximation. This relation can be intuitively explained as follows. The quantity U'(4>) 
changes due to the Brownian motion of in a manner that is independent of E, because the 
dependence of U'(4>) on E appears only at a higher order of E~ l . Equation ([36]) leads to an 
equation for the relaxation of the second element: 

(p 2 (t) 2 )^(p 2 (0f)-2 1 - 1 U^t. (37) 

It is noted that this form is independent of the temperature of the heat bath. 

Fluctuations around the average relaxation path, which depend on the temperature, can be 
straightforwardly estimated from ( |3~2"D as follows: 



((E - (E)) 2 ) ~ 7 " 2 /* t &t'&t"{U'{t')U'{t")){i{t')i{t")) ~ 2 1 - x TU' 2 t. (38) 
Jo Jo 

Thus, the distribution of E becomes broader with time, and this variance is proportional to 
the temperature T of the heat bath. The results represented by (|36"D and ([3g|) agree well with 
the numerical results given in §5. 

The ensemble average for the relaxation time t r is obtained from eq. ([36]) 



M - (39) 

This is again independent of the temperature. Thus the discussion given in §4.1 is also valid at 
finite temperature. 
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5 Numerical Experiments 



In this section, we check the validity of our analytical approximations, using results obtained 
by numerically integrating the Langevin equation (0) directly. In order to demonstrate the 
generality of our results, we adopted three forms of the potential U with < U < 1 (\5U\ ~ 1): 

i(l-cos0), 

^{2 - (cos# + cos g9)}, 

1 1 

cosh(acos^) cosh(a) 

Here, g is the golden mean (1 + \/5)/2, and the parameter a represents the steepness of the 
potential U$(6). U\{9) and U 3 {9) are periodic functions with period 2n, while U2(9) is quasi- 
periodic with two incommensurate frequencies. U 3 (9) possesses a steeper form as a increases, 
and it approaches a train of (^-functions with period 2tt as a — > oo. We set a = 10 to realize 
a sufficiently steep shape of the potential. The values of 5U 2 and U' 2 are given in Table I for 
each potential form. 

The relaxation at a finite temperature T was studied by taking an ensemble average. We 
integrated @ for almost 1000 samples with random seeds for each set of parameters, and 
thereby obtained a distribution pt{x) of the variable x(t), and thus an ensemble average (x(t)) 
over this distribution. For every sample in a given ensemble, we used the same initial conditions 
(given below), and thus the differences among paths for individual samples resulted only from 
the effect of the random force. 

In this section, we study the scaling form for the relaxation of the total energy E (not p 2 ) 
by considering variations of 7, U (6) and T. For this purpose, we chose an initial conditions at 
a given value of T with initial energy Eq as 

#i(0) = 0a(O) = pi(0) = 0, p 2 (0) = ^2E - T. (40) 

Here, Eq is a constant value, independent of T. By using these initial conditions, the total energy 
E was evaluated with the initial energy Eq for a number of temperatures, after a time to allow 
for the relaxation of the first element. After this relaxation, it is obtained that p?/2 ~ T/2 and 
Pl/2 ~Eq-T/2-U (i.e., E = p 2 /2 + p 2 2 /2 + U~E ). 

5.1 Underdamped case 

We first give the numerical results for T = 0. Figure [l] displays the energy relaxation for the 
three kinds of potentials described in Table I. For the underdamped case, the values 7„ in Table 
I were used for each potential. The relaxation profile obtained from the numerical experiments 
is fit well by E 2 = Eq — 7 5U 2 t, for 5U 2 corresponding to each potential form. Thus the 
numerical results confirm the analytical expression fl20P . 

As a next step, we examine the relaxation process at finite temperature. In Fig.|2|(a), 
ensemble averages (E 2 ) are displayed as functions of time t during the relaxation process. 
This figure includes numerical results from six different set-ups, using the two potential forms 
U\ and U3 each with three different temperature T = 0.01, 0.1 and 0.5. For a small time t, the 



U x {9) = 
U 2 (9) = 
U 3 (9) = 
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means (E 2 ) exhibit the same scaling relation as in the T = case: (E 2 ) ~ Eq — 7 5U 2 1. This 
is consistent with (0). As the time increases, the mean begins to deviate slightly from the line 
defined by (|25|). At higher temperature, there is a tendency to faster relaxation. 

In these computations, the variance of E 2 (i.e. ((P 2 — (E 2 )) 2 )) was found to increase linearly 
with time, as is seen in Fig.[|(b). This figure includes the variances for the above mentioned six 
cases, as well as those for two more cases with a different value of E Q at T = 0.01 for U\ and U3. 
The last two sets were included to examine the effect of changing the energy E . From these 
results, the variance is shown to satisfy the scaling relation ((E 2 — (E 2 )) 2 )/E T ~ A'~f5U 2 t, as 
expected from the analytical result in §5 (|27|). Note the slight deviation upward at large time. 



5.2 Overdamped case 

In this subsection, we report the numerical results for the relaxational behavior in the over- 
damped case. Here the parameter values 7 G in Table I were used. When T = 0, we found 
that E(t) beginning from a high-energy state relaxes as in Fig|| By considering a variety of 
potential forms (see Fig||), the result E ~ E — 7 _1 U' 2 1 is obtained, confirming the analytical 
result, (|33|). 

For finite temperature T, the ensemble average was studied numerically, and the same 
relaxation form, (E) ~ E — 7 _1 U' 2 1, is found to hold within the range of temperatures 
considered, as displayed in Fig.|](a). This figure includes results for six parameter sets: U\ and 
U 3 each with T = 0.01, 0.1, 0.5. 

The variance of an ensemble exhibits a linear increase with time, as seen in Fig.[|(b). From 
the numerical results, it obeys the scaling form ((E — (E)) 2 )/T ~ 2 7 _1 U' 2 1, which agrees with 
the analytical relation fl38|). We see that there is slightly larger deviation from the expected 
scaling relation here compared to the small deviation exhibited by the mean value in Fig.|](a). 



6 Crossover Regime 

In this section we study the dissipation of energy in intermediate values of damping between the 
overdamped and underdamped limits. In the preceding sections, the relaxation of the ensemble 
average of E (or E 2 ) was shown to be independent of the temperature. For this reason, we 
confine our study of the crossover regime to the case T = 0. The result we obtain here is 
expected to hold at a finite temperature also. 

Recall the relaxational behavior at T = obtained for the underdamped and overdamped 

limit: 

f 7 #7* 1 

AE 1 -v*- 1 . (41) 

_ 7 -lp> ; p < 7. 



dt 



Here E ~ P 2 /2 and U <C E. In Appendix 2, we obtain more general approximations of p± 
and E, for any values of 7 and P. Here we give an intuitive description of the change in the 
relaxational behavior with the decrease of P for a fixed value of 7. 

The first element is driven by the second element and damped by the frictional effect of the 
heat bath. When 7 is sufficiently small, and P is sufficiently larger than 0(7), this damping 
due to friction is negligible, and the dynamics of the first element are approximately determined 
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solely by its interaction with the second element. This leads to p\ ~ SU/P. From this we obtain 
the upper form in (f41"l). As the momentum P decreases, the interaction between the first and 
the second element becomes stronger, and this is accompanied by the growth of p%. Due to this 
input of energy to the first element, it rapidly dissipates energy to the heat bath. In such a 
situation, the value of p\ cannot be estimated as SU / P, due to the large friction. Instead, it is 
estimated as 7 _1 t/' for sufficiently small values of P. The relaxation equation in this case thus 
assumes the lower form in (|4lD . 

From the above picture, we expect a smooth interpolation between the two relaxation forms 
in fl4ip . The effect of the heat bath on the first element is comparable to that of the second 
element when the energy E (~ P 2 /2) satisfies the condition >y5U 2 E- 1 /2 ~ 7 1 £7 72 . Taking this 
into consideration, we define the crossover energy for the relaxation behavior, 



7 2 5U 2 

E c = -U^. (42) 
2U' 2 

This represents the energy at which the two relaxational forms are equal. By introducing t c 
defined by 

tc = I =r, 43 
2U> 2 

the normalized energy e = E/E c , and the normalized time r = t/t c , ([O) can be written as 

-2e"\ £»1, 



d£ 

dr 



(44) 

£ < 1. 



These equations include neither 7 nor U explicitly. Hence, the value of de/dr is expected to 
depend only on the value of e, regardless of the value of 7 and the shape of U. 

The relaxation form e obtained numerically is displayed in Fig.|5| for the case T = 0. To 
obtain this figure, we computed (0) at T = for various values of 7 and plotted the relation 
between e and de/dr after the first element reached thermal equilibrium. We see that the 
value of de/dr becomes proportional to e^ 1 as e increases beyond 1, while de/dr approaches 
— 1 and becomes independent of the value of e for smaller e, as expressed by (|44"1). This figure 
shows that the relaxation form changes smoothly between the two limiting regimes. We also 
remark that the nature of potential form U (8) seems to affect the sharpness of the crossover 
near E = E r . 



7 Summary and Discussion 

In this paper, we have studied the relaxational behavior from a far-from-equilibrium (highly 
excited) state for a system consisting of internal and external elements in which only the external 
element is in contact with a heat bath. In spite of the simplicity of the model, our system 
exhibits a power-law-type slow relaxation when the potential is bounded, i.e., for \U\ < 00. 
This slow relaxation is observed when the inner part of the system is excited to a high-energy 
state, with energy much greater than \8U\. The slow relaxation of the internal part results from 
its indirect contact with the heat bath. In spite of the system's contact with the heat bath, 
guaranteeing its approach to equilibrium, the observed slow relaxation is in a strong contrast 
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with the rapid exponential decay exhibited by systems in full contact with a heat bath. We 
therefore find that a Hamiltonian system in partial contact with a heat bath can maintain its 
own internal dynamics for a much longer time than one in full contact with a heat bath. If a 
system has an internal part further separated from the heat bath, the Hamiltonian dynamics 
of this part can be preserved much longer. 

Suppose the momentum of the second element P2 is constant in time. Then the present 
model becomes + 70 + 2£7'(0) = — 7£>2 + The Brownian motion of a particle in a periodic 
potential U((f>) in thermal equilibrium has been analyzed using the Fokker-Planck equation for 



this model (fLO, [TT[ in the contexts of Josephson tunneling junction [12 1, super-ionic conductors 



[PL3| , and two-mode lasers |14j| . In contrast, we are interested in the autonomous relaxation of 
P2 from a highly excited state, through the internal dynamics obeying (0). This autonomous 
change of P2, due to the internal dynamics of the composite system, brings about slow relaxation 
characterized by power laws, even though the system is neither subject to special conditions, 
such as being placed at a phase transition point, nor exists in a glassy state involving many 
degrees of freedom. 

While our interest in Hamiltonian systems in partial contact with a heat bath originates in 
our desire to gain an understanding of the long-term maintenance of high-energy states, such 



systems have also been studied in relationship with thermal conduction ||15|| , for example, in the 
situation that two ends of a chain are in contact with two heat baths at different temperatures. 
According to several simulations performed on systems of this type, the Fourier law of thermal 
conduction is not always obeyed. It appears that certain conditions concerning the internal 
dynamics, as well as a suitable relation of the system with the heat baths, are required in 
addition to ergodicity (or mixing), in order for the Fourier law to be obeyed. These studies also 
suggest that parts of systems in indirect contact with heat baths can maintain over extended 
periods of time the characteristic dynamics of isolated Hamiltonian systems. 

In a class of Hamiltonian systems in the absence of a heat bath, very slow relaxation 
to equilibrium was suggested in the form of the 'Boltzmann- Jeans conjecture' 0, and such 
behavior was reported in Refs. H, 0|. Long-term storage of a large amount of energy within 
a degree of freedom through this mechanism has also been studied in a Hamiltonian system 
with chaotic behavior ||, in which localized energy is transferred to other degrees of freedom 
only after a very long time. In that case, the energy dissipates logarithmically in time. We 
have reported in this paper that situations exist in which localized energy in a system dissipates 
slowly, according to power laws, and that long-term energy storage is possible even if the system 
is integrable and in contact with a heat bath. From this point of view, our result can be regarded 
as an extension of the Boltzmann- Jeans conjecture to include dissipation, while the relaxational 
behavior itself is modified from the logarithmic form of the Boltzmann- Jeans conjecture. As 
mentioned in the Introduction, long-term energy storage has been reported in some protein 
molecules 0, and the relevance of this storage to molecular motors (for muscle contraction) 
has been discussed. To apply the present ideas regarding slow relaxation to proteins, we need 
to examine relaxation from a high-energy states for a system with more degrees of freedom. 
This problem will be studied in a future publication |I6| . The main question to be addressed 
is whether the relaxational behavior exhibits a power-law-form or a logarithmic form when 
Hamiltonian systems more complex than that considered in the present study are in partial 
contact with a heat bath. 

This work is supported by Grants-in-Aid for Scientific Research from the Ministry of Edu- 
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Research Resources of Japan Atomic Energy Research Institute. 



8 Appendix 1 

In this appendix, an expression for the ensemble average (pf) is derived from fll7|) . Since this 
equation is a Langevin equation with multiplicative noise, the drift term of the corresponding 
Fokker-Planck equation is not necessarily — / y5U 2 . We carry out the Kramers-Moyal expansion 
[|Tofl for the evolution of the variable z = p\ to obtain the Fokker-Planck equation for the 
probability P{z,t) in the form dP(z,t)/dt = dD^P/dz + (l/2)d / dzD^ d / dzP . Computing 
straightforwardly from fll7|), we obtain the drift term 



and the diffusion term 



V2pl Pi 



P2 



where the second term of comes from the multiplicative noise. 

Since we confine our study to the situation with sufficiently large pi and sufficiently small 
T, this second term of can be ignored. Thus we have = —"fSU 2 to the order of 
approximation of the current analysis. Hence, the ensemble average (p^SU^) can be decomposed 
as (p- 2 SU)(0- 



9 Appendix 2 



In this appendix, we give another derivation of the relaxation equation and a more precise 
expression of the crossover condition at T = 0. Solving the equation of motion for the first 
element ( p\ = —dU(<f>(t))/d<f> — jpi), we obtain 



pi(t)= r dT d£(^i)) e -,, 

We write the periodic potential U ((f)) as a Fourier series: 

' u = H a i c os(^i(0-0i)) ; 
j 

U' = -^a^sin^^-^-)), 



(45) 



(46) 



where we assume U = without loss of generality. Then, assuming that cj)(t) = <pt, (f45|) can 
be written, 



Pi = ~ 1^ a j"j / dr sm{uj<j){t - r) + <^)e 7 = ^ 



(47) 
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Here Uj and f/j are the jth components of the expressions in ([R]). This solution implies that 
Pi = (p^U, E = -7 4>~ 2 TP (<j>Uj > 7 for all j), 

and 

Pi = l^U, E = -7" 1 I7 7 2 ^ ( j )UJ . ^ 7 f or a u 

These equations are equivalent to those obtained in §6, but in the presents derivation, the 
crossover condition is more explicitly represented as <pujj ~ 7. Since Uf = oo 2 -U' 2 , this crossover 
condition can be expressed as 

d> ~ 7 =^=, 

which is consistent with the condition discussed in 56. 
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Table 1: Characteristic values 5U 2 and U' 2 for the three potentials. 7„ and 7 are the values of 
the relaxation coefficient 7 used in the numerical computations for the underdamped and the 
overdamped limit. 



Figure 1: Energy dissipation for T = 0. Relaxation from the initial conditions ([40]) with E = 10 
were computed. E 2 is plotted versus the scaled time 7 5U 2 t. The values of 7 are those given 
for 7„ in Table I for each potential. With the scaled time 7 5U 2 1, the relaxation curves for the 
three different potentials coincide. This is consistent with the analytical estimate. 



Figure 2: (a) Relaxation process of the ensemble average (E 2 ) for T = 0.01, 0.1 and 0.5, for each 
of the potentials U\ and U3. The values of 7 are those of 7„ given in Table I, while the initial 
conditions are given by ( [40]) with E = 10. The scaled time r y5U 2 t is used for each potential. 
The line indicates the analytical approximation represented by (^). (b) Temporal evolution 
of the variance of E 2 corresponding to the six cases in (a) and also for two additional cases 
with Eq = 12 and T = 0.01 for U\ and U3. The line indicates the analytical approximation 
represented by (pf). 



Figure 3: Energy dissipation for T = 0. The value of 7 fpr each potential is the corresponding 
value of 7o given in Table I, and the initial conditions are given by (|40|) with E = 10. With 
the scaled time 7 -1 U' 2 1, the relaxation curves for the three different potential coincide, in 
agreement with the analytical approximation. 



Figure 4: (a) Relaxation process of the ensemble average (E) for T = 0.01, 0.1 and 0.5, starting 
from the initial conditions ( |40"D with E = 10. The potential forms used are U\ and U$, while 
the parameter values 7 are given by 7 Q in Table I, corresponding to the overdamped case. 
The scaled time 7 _1 U' 2 1 is for each potential. The line indicates the analytical approximation 
represented by (36). (b) Temporal evolution of the variance of E corresponding to the six cases 
in (a). The line indicates the analytical approximation represented by (|38|). 



Figure 5: Dependence of E on E for a variety of values of 7, which range from the underdamped 
to overdamped case. Here E c and t c are determined for each potential according to the value 
of 7 (see (f|2]) and (^)). The functional forms of E(E) for the two potentials Ui and U3 are 
essentially identical, except in the intermediate range around E/E c = 1. 
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